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Abstract 

We study the Einstein-Klein-Gordon equations for a convex positive 
potential in a Bianchi I, a Bianchi III and a Kantowski-Sachs universe. 
After analysing the inherent properties of the system of differential 
equations, the study of the asymptotic behaviors of the solutions and 
their stability is done for an exponential potential. The results are 
compared with those of Burd and Barrow. In contrast with their re- 
sults, we show that for the BI case isotropy can be reached without 
inflation and we find new critical points which lead to new exact so- 
lutions. On the other hand we recover the result of Burd and Barrow 
that if inflation occurs then isotropy is always reached. The numer- 
ical integration is also done and all the asymptotical behaviors are 
confirmed. 

PACS number(s): 98.80.Cq, 98.80.Hw 
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1 Introduction 



Inflation, first introduced by Guth [|IJ, was introduced in the standard cosmo- 
logical model to solve the homogeneity, the isotropy and the horizon problem. 
The latter is well explained due to the fact that inflation is characterized by 
an exponential or power-law expansion of the universe and at the same time 
a quasi-constant behavior of the Hubble horizon. 

On the other hand the homogeneity and the isotropy problem is in fact 
not well explained because from the start the homogeneous and isotropic 
Friedman-Lemaitre metric is used. To really solve the problem one should 
start with an arbitrary metric and show that inflation takes place and that 
the universe evolves towards a Friedman-Lemaitre metric. The problem of 
the onset of inflation was considered numerically for spherical inhomogeneous 
cosmologies |2|, |3[] and a semi-numerical analysis was done for inhomogeneous, 
quasi-isotropic universes H using the long wavelength iteration scheme || |6| . 
They showed that a large initial inhomogeneity supresses the inflation stage. 
Because of the analytical difficulties of the task, one can as a first step, 
consider only a homogeneous but anisotropic metric and try to solve the 
isotropy problem. 

This task was first undertaken by Collins and Hawking pj who showed 
that, within the Bianchi type universe filled by matter satisfying the domi- 
nant energy condition and positive pressure criterion, the isotropy problem 
can only be solved for the types I, V, VII Q and VII^. They showed also that 
only a subclass of vanishing measure in the space of all homogeneous initial 
conditions can approach isotropy. 

With the presence of an inflationary stage, when the dominant energy 
condition is violated, there was hope to obtain a cosmic no-hair theorem. 
The study done by Heusler || showed there is no no-hair theorem for a 
real scalar field having a convex positive potential with a vanishing local 
minimum in a Bianchi type universe. In fact isotropy can only be approached 
if the underlying Lie group of the Bianchi type metric is compatible with a 
Friedman-Lemaitre model. 

Beside the Bianchi type metrics the Kantowski-Sachs model also describes 
a spatially homogeneous universe. This model with a perfect fluid description 
of matter and with or without a cosmological constant has been studied by a 
number of authors [[], [IT], [IT], [T^, [Tj|. They found an anisotropic asymptotical 
behavior of the model. Burd and Barrow analyzed this system with a 
real scalar field having an exponential potential as a source. Among others 
they found an anisotropic asymptotical behavior or that inflation can occur 
depending on the value of the coupling constant entering in the definition 
of the potential. The exponential potential is motivated by the fact that it 
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can be obtained, for example, by dimensional reduction of more fundamental 
theories [15| or in conformal equivalent theories of gravity whose Lagrangian 
is an arbitrary analytic function of the scalar curvature TO , |TT| , |TE |. 

In this work, we will consider a real scalar field with a convex positive 
potential, but not necessarily with a local minimum, in a Bianchi I, a Bianchi 
III and a Kantowski-Sachs model. These three model are in fact very closely 
related since it will be shown in section |2| that the differential equations 
describing the evolution of the models are the same and that only the con- 
straint equations differ. In section |3| we study the inherent properties of the 
differential equations and we recover in the case of a local minimum of the 
potential a result of Heusler || . We also show why the solution of the dynam- 
ical system must be studied for a potential having a vanishing potential at 
infinity. Section [| is devoted to the special case of an exponential potential. 
The result obtained will be compared with the analysis done in Ref . [frjj . In 
section |^ the numerical integration of the system is done and in section ^| the 
results will be summarized. 



2 Basic equations 

We consider a scalar field with a convex positive potential V in a homoge- 
neous universe having the following metric 

d s 2 = g ^e v = -e°e° + s^ew, (i) 

where 

e° = dt , e l = a(t) dr,e 2 = b(t) m , e 3 = b{t)f(#)d<f> , 

with /($) = f($) = sinhifi) or /(■#) = sini^d). We have respectively a 
Bianchi-I (B-I), a Bianchi-III (B-III), and a Kantowski-Sachs (KS) metric. 
The action is given by 

S = J (-J^q + \^e^ + v^9 d 4 x , (2) 

where g is the determinant of the metric and e M is the dual basis of 9^. By 
varying the action with respect to the metric we get the Einstein equations 

2H a H b + H 2 b +^ 2 = 8nG Q^ 2 + Vfo>)) , (3) 

2H b + 3H 2 b + A = 8 nG (-^ 2 + , (4) 

H a + H b + H 2 a + H a H b + H 2 b = 8nG (-^ 2 + V{ip)\ , (5) 
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with 



u h u b i ld2 f ret\ 

H * = a> H > = V k= —fW (6) 

For the B-I, the B-III and the KS cases, we get respectively k — 0, —1, 1. 

The Klein-Gordon equation is obtained by varying the action with respect 
to the scalar field. We get 

dV 

<p + {H a + 2H b )v + — = Q. (7) 

dip 

In the following we will express every quantity in units of the Planck mass 
m p = l/y&nG. This can be achieved by setting 8ttG = 1 in the Einstein 
equations. The system is fully determined by the independent equations 
(3), (4) and (|7|). One can easily show that equation (5) follows from the 
others. This is a reflection of the Bianchi identities. After some algebraic 
manipulations, these equations can be written as a set of four first order 
differential equations, which are k independent, and a constraint (eq.(3)) 
which is conserved in the evolution. 

This set of equations can also be written in function of the expansion rate 
9 and the shear tensor <j pu of the hypersurface of constant time S. We get 
for the first order equations 

e = - l -e 2 -2a 2 + v^)-^\ (8) 



a 



3^ V3 V3 

cp = iP, (10) 

* = ^+-% ' < U > 

where a = \ a^o^ = ^/T/3{H a - H h ) and 9 = H a + 2H b . ip is defined by 
eq. (10). As for the constraint equation, it reads 

^ 2 + p = ^ 2 + ^) + ^ 2 - (12) 

The solutions of eqs.(8)-(ll) do not depend on the type of the homogeneous 
model considered. Since eq.(|l~2D is conserved, we have to specify the homo- 
geneous model only in the initial conditions. In the four dimensional space 
with coordinates (9 , a , z = VV , ip) the B-I solutions are on the "lightcone" , 
the B-III ones are in it and the KS solutions remains outside the lightcone 
(see Fig. 1). 
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The advantage of taking these geometrical meaningful quantities is ob- 
vious: the shear tensor expresses the direction dependent deviation from 
the global expansion. Hence, a measures the anisotropy. As a criterium of 



isotropization we will not use the vanishing of a (as in |19|), but a stronger 
condition [0, Q which is 

a 

— — > , as t —>• oo . (13) 
9 

In this paper, we want to analyze the evolution of an expanding universe. 
Hence, we will restrict ourselves to positive values of 9. We want to know 
under which conditions the isotropization of this model is generic and when 
inflation occurs. In some cases, general properties of the solutions can be 
found by just analyzing the differential equations. 



3 Inherent properties of the dynamical sys- 
tem 

This section generalizes, in some respect, a result obtained by M. Heusler 
||. In the B-I (k = 0) and the B-III (k = 1) case, the asymptotic behavior 
of the solutions can be obtained directly by the study of eqs.(8)-(0). Since 
we are interested in the evolution of an expanding universe, we suppose that 
6 is positive or equal to zero at some time to. We also restrict ourselves to 
positive and convex function V. From eq.(8) and eq. (|T^) , we get the following 
inequality relation (k < 0) 

9<V-^9 2 <0. (14) 
Using the last inequality relation and again eqs.(8jT2|), it is easy to obtain 

e + 6 2 = 3V-^>0. (15) 
tr 

Setting V = ab 2 which is positive, we can write 

V 



By hypothesis, V is positive at t . Substituting eq.fllq) in eq. (|15|) , we have 
that V/V > 0. Therefore, V can only increase with time and 9 will remain 
positive after t . We have shown that 6 is a positive, monotonic decreasing 
function and thus, we must have 

and 9 -> 6^ > as t -> oo . (17) 
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This function converges towards the critical value 9^. It is also easy to get 

^ = p-3a 2 -^ 2 <0. (18) 

Since the l.h.s of the inequality is a sum of negative terms, we have that each 
of them must vanish as t — ■> oo. Using this fact in the constraint equation, 
we obtain 

9 -> VW as t -> oo . (19) 

Eq.(ll) describes a damped harmonic oscillator, so we know that the system 
must "land" at the minimum of the potential Vq and thus 



#oc = pV . (20) 

For an exponential potential, the only minimum is at infinity and is value is 
zero. Let us summarize the results by the following theorem 

Theorem 3.1: Let befc = 0orfc = — 1 (Bianchi-I and Bianchi-III case). 
If at a given time to we have an expanding universe, that is 9(to) > 0, then 

(i) 9(t) < 0, 9{t) > for all t > t , 

(ii) a(t), if)(t), — > 0, for t — ► oo and 



(iii) 9 — > yoK)) where Vq is the minimum of the potential and in particular 
9^0 for an exponentially potential (i.e. V = Voe~ Xlf ) for t — > oo. 

This theorem restates the proposition 1 in |§ except than it can now also be 
applied for a convex potential V having its minimum at infinity. If Vq > 0, we 
can use the known result that the B-I and the B-III converge exponentially 
to the isotropic De-Sitter universe ||20|| . If Vq = at some finite value ifo then 



we only have isotropization in the B-I case ||. 

The above theorem gives no answer on isotropization for an exponential 
potential in a B-I or a B-III model. Indeed, 9 and a vanish together. In this 
case, the asymptotic solutions of the dynamical system (eqs(8)-(ll)) must be 
found to know under which conditions isotropy and/or inflation can occur. 
Since the equations of the dynamical system do not depend on k, the KS case 
will also be solved. This system has been actually studied in a paper of Burd 
and Barrow |14|]. However, their analysis on the system was not complete. 



Their analysis covers a part of the phase space only. Indeed, only two of the 
critical points were found and some conclusions turn out to be incorrect. 

In the next section, the analysis of the exponential potential case will be 
done in detail and we will compare our results with Ref. [JTJj]. 
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4 Study of the dynamical system with an ex- 
ponential potential 

4.1 The singular points 

For a qualitative discussion of the system of differential equations (8)-(ll) we 
have to determine the asymptotic behavior near the critical points. The usual 
methods of treating the problem rapidly turn out to be inadequate for this 
case: the singular points are highly non-hyperbolic so that the linearization 
or even the transformation to a normal-form lead nowhere. Usually, the 
constraint equation can be used as a Ljapunov-function. But in this case, 
eq. flT2| ) has no isolated root and therefore is of no help. There is no other 
obvious candidates for a Ljapunov-function. 

The problems originate from the fact that there are only quadratic terms 
on the right hand side of eqs.(8)-(ll). There is a standard procedure for 
such systems: one has to "divide" all variables but one by the remaining 
variable and rewrite the system in terms of these new variables. Of course 
this transformation is only non-singular if the divisor does not vanish in the 
whole region of definition. In Ref. ||14j| , the remaining variable was chosen to 
bey = V($), which does not vanish in finite regions of the phase space. Here, 
we prefer to distinguish the variable 6 since it will allow us to find additional 
critical points (in finite regions of phase space). This transformation is well 
defined outside the origin. 

Defining the variables S, U and P by 



a = se, z = ue, v = po 



(21) 



and the time r by 



d 

dr 



1 d 
6di 



(22) 



the differential equations (8)- (11) transform into 



& = 9{-- - 2S 2 + U 2 - P 2 ) 



(23) 



S' 



1 




P 2 + 2S 3 - SU 2 + SP 2 



3V3 



(24) 



P' 




(25) 



(26) 
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and the constraint equation is 



S2 + lj2 + \ p2 + ^ = \- (27) 

The dynamical system defined by eqs.(24)-(26), which are independent of 9, 
will be referenced from now on by the tag (*). 

For the Bianchi I case (k = 0), equation (27) describes the surface of 
an ellipsoid, say E, which separates the Kantowski-Sachs (inside) from the 
Bianchi III (outside) solutions. The Bianchi I model (k = 0) has to be treated 
carefully because in this case the system can only evolve on a two dimensional 
submanifold W in the space of the variables (S, U, P), defined by 

W = f-\{ 1 -}) with f(S,U,P) = S 2 + U 2 + ~P 2 . (28) 

It may happen that W coincides with the stable manifold of a critical point 
which is unstable with respect to the whole system (*). In this case the point 
is nevertheless stable for the Bianchi I metric. 

We will concentrate on determining the asymptotic behavior of this re- 
duced system. 9 is then given by a simple integration of (23). After a short 
calculation we find for the system (*) three discrete critical points and a 
one-parameter family of critical points: 

Pi : S = U = 0,P = (29) 

P = - for A < V6 (30) 
^ P = ^ (3D 



P 2 : 


s = 


V 18 ' 




s = 


1 2-A* 
2v/3A 2 + l' 


£ : 


u = 


0, 3S 2 + Ip 2 = . 



x/2(A 2 + l)' A 2 + l 
1 with S E y/l/3] ■ (32) 

Comparing these critical points with the results of Burd and Barrow [[HJ] one 
easily verifies that P<i and P3 correspond to their critical points (IV) and 
(III) whereas Pi and the points of £ have no counterpart in their paper. 
Using the variables in Ref. pi, Pi and £ can be found as critical points 



lying at infinity in their phase space. 

To be sure that there are no additional critical points at infinity in our 
coordinate system either, we first introduce spherical coordinates 

S = r sin??, U = r cos ip cos P = r simp cos'd (33) 
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and then map the points at infinity (r = oo) to the surface of a unit sphere 
by the transformations 

r — ► p = dr — > dr) = . (34) 

r + 1 1 — p 



At infinity (p = 1) we thus find 

2-cos 2 tf(l + 2cosV) (35) 



dr] 

^ cos <£>cos$(2 cos 2 cp + sin 2 cp) (36) 



dr] 2 
d$ cos 

— = — 1= (2 — cos 2 ?? + cos 2 cp cos 2 r3 — V%\ sin <p cos 2 o? sin cos r3) . 

dr] 2^/3 ^ 

For a critical point at infinity, expressions (35) - (37) have to vanish simul- 
taneously which is easily seen to be impossible. Therefore the points Pi, P2, 
P3 and the points of £ are really the only critical points of the system (*). 

It is worthwhile to note that the critical points correspond to exact solu- 
tions of our original system ((8)- (11)). Indeed, if we consider the pull-back 
of a critical point (Sk, Uk, Pk) of (*) for a non-constant potential (A > 0), the 
volume expansion 9 is given by 

9 = 9 e AkT (38) 

with A k = —1/3 — 2S 2 + U 2 — P 2 \(s,u,p)=(s k ,u k ,p k )i an d as a function of the 
original time t we obtain the following exact solutions 

-h— h *= ~h vm =w • (39) 

The stability analysis of the critical points will give the stability of the above 
exact solutions. 



4.2 Stability analysis of the singular points 
4.2.1 The Critical Point Pi 

The linearization of (*) around the critical point P\ : (S, U, P) = (— (2\f?>)~ 1 , 
0,0) is diagonal and has eigenvalues £g,p = —1/2 and Ey = 1/2. So Pi is 
unstable both to the future and to the past. Starting around Pi, the critical 
point can never be reached by any solutions of the dynamical system. 

Using eq. (|39|) , the critical point Pi will correspond to the unstable 
solution 

= \>°= "7y=Y i> = °> ^( $ ) = ■ ( 4 °) 
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4.2.2 The Critical Point P 2 

In P 2 the linearization of (*) has the eigenvalues E\ — — 1 + A 2 /6 (twice) and 
62 = — 2/3 + A 2 /3, so the critical point is asymptotically stable (with respect 
to all the system (*) and for r — > 00) as long as A < y/2 and unstable for 
y/2 < A < y/E. 

Since P2 lies on the ellipsoid which separates the KS from the Bill do- 
main, its neighborhood intersects all three types of universe. More precisely, 
starting from a point P, near P2, defined by 



S = 5S, U=^-^ + 5U, P = ^ + SP (41) 

the constraints equation (27) expanded to first order tells us to which universe 
P belongs since we have 



Depending from which point we start the integration, the solution can be in 
either type of universe. 

As mentioned earlier, the Bianchi I (k = 0) metric has to be consid- 
ered very carefully because the system is constrained by (27) to the two- 
dimensional submanifold W. The tangential space to W in the point P2 is 
given by 

T P2 W = {X G R 3 : X JL V/(P 2 )} • (43) 



But V/(P 2 ) = (0, y 2(6 — A 2 )/3, A/3) is an eigenvector for the eigenvalue £2, 
and so Tp 2 W coincides with the eigenspace to the eigenvalue e±. That is, W 
is the stable submanifold through P 2 for A < y/E and consequently P 2 is in 
all this range an asymptotically stable critical point for Bianchi I solutions, 
contrary to the Bianchi III and Kantowski-Sachs solutions. This particular 
behavior for the BI case was overseen in Ref . ||14|| . 

For the special values A = a/2 and A = y/E at least one eigenvalue vanishes 
and the linearization does not contain enough information to determine the 
stability of P 2 . By projecting the system (*) on a center manifold through 
P 2 |2l], Q we find that P 2 is unstable for A = y/2. In the case A = y/E 
the center manifold is two dimensional and therefore the stability analysis 
is much more difficult. We did not found the asymptotical properties of the 
solutions around this singular point for this value of A. 

The critical point P 2 transforms back to the exact solution 

6 2 , x 2(6 -A 2 ) 
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For A < \/2 this solution violates the dominant energy condition (tp 2 — V < 0) 
and thus inflation occurs. This can be directly seen since from the above 
equation we obtain 

a = a t 2/x2 (45) 

For the asymptotically stable range (A < y/2 for KS, Bill and A < a/6 for 
BI) of P 2 the solution is also asymptotically stable in the sense that solutions 
which start near it converge to it. 



4.2.3 The Critical Point P 3 

The eigenvalues of the linearization of (*) in P3 are 



A 2 + 2 -(A 2 + 2)±^(18-7A 2 )(A 2 + 2) 

£l ~ ~2(A 2 + 1)' £2 ' 3 ~ W^T) ' ( } 

For A < V% all eigenvalues are real, e 2 and e 3 have different signs, so P 3 is 
unstable. For A = a/2, P 3 coincides with P 2 so we conclude as before that 
the critical point is unstable. In the range a/2 < A < J 18/7 all eigenvalues 



are real and negative, while for A > y 18/7 e 2 and e 3 become complex with 

the same negative real part. So P3 is asymptotically stable for A > v2. 
Again we find the following exact solution after the pull-back of P3 

which is asymptotically stable for A > a/2. 



4.2.4 The Critical Points of £ 

Since the critical points in S are on the ellipsoid E, we find as before (section 
4.2 .2|) that there are points belonging to all three type of universe near S. 

For the linearization of (*) around a point Ps = So, 0, Pq of S we find the 
eigenvalues 

ei = 0, e 2 = 1-^P , e 3 = ^(2 + V3So) • (48) 

If P > 2/A, e 2 and e 3 have different signs and the critical point is unstable 
both to the future and to the past. 

If Pq < 2/A, e 2 and e 3 are both positive, so Ps is unstable to the future. 
The stability for r — > — 00 is determined by the behavior of the system on a 
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center manifold |2I, ^2|. It is easy to see that S itself is a center manifold 
and that the restriction of (*) to S is just the trivial system 



S' = 0, U' = 0, P' = . (49) 

That means that for r — > — oo P^ is an attractor. It his also easy to see 
that for A < a/6 E as a whole is a past-attractor for (*). Indeed, £2 is then 
positive for all points Ps- 

The pull-back of points out of S gives an exact solution of the form 

6> = i a = ^ = ^, V(S) = 0, (50) 
with the condition | — 2 <Tq — ^ = 0- 



4.3 Summary 

The cases with A = (constant potential) can be treated analogously. For 



the critical points P2 and P3 this has been done by Burd and Barrow \TA 
while for Pi and S the results founded are still valid since the asymptotic 
behavior does not depend on A. 

The results of this section are summarized in table 1 (isolated critical 
points) and table 2 (family E). They show for every critical point the type of 
the corresponding exact solution (BI /III: Bianchi I /III or K-S: Kantowski- 
Sachs), its Isotropization (Is) and inflation (In) (e: exponential inflation, p: 
power-law inflation or -: no inflation), and the stability (St) of the critical 
points (|: asymptotic future-stable, t: asymptotic future-stable for Bianchi I 
solutions and unstable otherwise, Is' past-convergence to £ and -: unstable). 

From table 1, it follows that for A < a/2 the unique asymptotical stable 
solution is represented by the point P2. This solution describes an inflationary 
and isotropic universe. 

For a/2 < A < a/6, isotropy can be reached without inflation if we restrict 
ourselves to the BI universe. This fact is compatible with a Collins and 
Hawking result 0, which states that for ordinary matter, within the BI 
universe, isotropy can be reached without inflation. Indeed, for a/2 < A < a/6 
the equation of state of the matter field corresponding to the P2 solutions is 
given by 

p = up (51) 

with uj g] — 1/3, 1[ and contains the ordinary matter case (0 < uo < 1). 

For A > a/6, the unique asymptotical stable solution is represented by 
the point P 3 . This Bill solution describes an anisotropic universe. 
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There is also a unique attractor in the past: the £ manifold. As a con- 
sequence, we must start the numerical integration near £ for having all the 
history of the evolution of the model. 

The next section will be devoted to the numerical integration. 



5 Numerical Results 

In this section we will follow numerically the evolution of the system (8)- 
(11). We have seen in section |] that the points of the surface £ are the only 
past-attractor of the system, so it seems obvious to study the development 
of solutions starting near the corresponding exact solution 

= \, <r = v j, V=^, V{$) = (52) 

with t — *■ and with 3<7q + 3/2?/>q = 1. Since these exact solutions are in the 
BI universe, the numerical solutions starting near £ can be in any three type 
of universe (e.g. BI, Bill or KS). 

To distinguish the different future-attractors it is convenient not to plot 
the variables 9 and a themselves (they always vanish as 1/t) but the quan- 
tities 9 ■ t and a ■ t which converge to constant values at the critical points. 

5.1 Kantowski-Sachs Solutions 

For the initial asymptotic behavior 

0(*o) = — , 0"(*o) = 7 , Uyto) = — , Wo) = T , 

H) H) H) K> 

we will have a Kantowski-Sachs solutions if the following inequality is satisfied 

2a Sa + ipoSip > . (54) 

Fig. 2 shows the results for some values of A less then the critical value \/2. 
The solution clearly converges towards the exact solution corresponding to 
the critical point Pi- The corresponding phase portraits are given in fig. 3. 

The development of a solution with A > \pl is shown in fig 4. The 
solution approaches again the (now unstable) attractor P 2 , is repelled and 
finally converges towards the exact solution corresponding to the attractor 
P 3 . 
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5.2 Bianchi III Solutions 



For a Bianchi III solution of the form ( |53"D the integration constants have to 
satisfy the condition 

2a 5a + i/) 6ip < . (55) 

The solutions with A < \[2 (fig. 5) initially tend to the Bianchi III solution 
corresponding to P\. But since this is an unstable attractor the system 
finally evolves towards the solution P 2 . Fig. 6 shows the corresponding phase 
portraits for these cases. 

For A > \/2 the system again approaches the unstable Pi-solution in the 
beginning, is repelled and converges towards the P3-solution which is the 
only attractor in this range of A- values (fig. 7, phase portraits: fig. 8). 

5.3 Bianchi I Solutions 

For a Bianchi I solution, the perturbations of the initial state (^) have to 
satisfy the equation 

2a 5a + ip 5?p = . (56) 

For Bianchi I solutions, the exact solution corresponding to P 2 is attrac- 
tive not only for A < \[2 but for the whole range for which P 2 is defined, that 
is for A < \/6 (fig. 9, phase portraits: fig. 10). The numerical calculation 
confirms our analysis for v2 < A < y/E: the solutions still isotropize but 
are no longer inflationary as one can see from fig. 11 where we have plotted 
ip 2 — V(<&) which violates the dominant energy condition when negative. 

For A > \/6 we have again the solution to P3 as the only attractor. An 
example for this case is shown in fig. 12. The system quickly evolves towards 
the now unstable solution P 2 , is repelled and converges towards the solution 
to P 3 . 

6 Conclusions 

We have studied the Einstein-Klein-Gordon (EKG) equations for a convex 
positive potential in a Bianchi I, Bianchi III and a Kantowsky-Sachs universe. 

After analyzing the inherent properties of the equations, it was shown in 
section [| why a detailed analysis of the solutions of the EKG equations was 
needed for a vanishing potential at infinity. By taking an exponential poten- 
tial (V = Voe~ Xv ) it was shown for which values of A inflation and/or isotropy 
where reached asymptotically. We recovered the results of Ref. jI3J but also 
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found new asymptotical behaviors and new exact solutions represented by 
the singular points Pi and the sub manifold S. 

We also found that for some values of A isotropy can be reached without 
inflation in a Bianchi I universe. But when inflation occurs (A < y/2) then 
isotropy is always reached. We also integrated the equations numerically 
to obtain the all history of the evolution of the model. All the founded 
asymptotical behavior were confirmed numerically. 
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Caption 



Table 1: 

Summary of results for isolated critical points 
Table 2: 

Summary of results for critical points in E 
Figure 1: 

Region of the space with coordinates (8, a, z, ■0) where the solutions lie as 
function of k. The solutions for the B-I case are restricted in a submanifold 
of codimension 1, the "lightcone". 

Figure 2: 

Evolution of Kantowski-Sachs solutions with A < \/2. 
Figure 3: 

The corresponding phase portraits for the KS case with A < \/2. 
Figure 4: 

Evolution of a Kantowski-Sachs solution with A > \/2. 
Figure 5: 

Evolution of Bianchi III solutions with A < \/2. 
Figure 6: 

The corresponding phase portraits for the Bill case with A < y/2. 
Figure 7: 

Evolution of Bianchi III solutions with A > y/2. 
Figure 8: 

The corresponding phase portraits for the Bill case with A > \/2. 
Figure 9: 

Evolution of Bianchi I solutions with A < \f&. 
Figure 10: 

The corresponding phase portraits for the BI case with A < a/6. 
Figure 11: 
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Inflation can only occur when the dominant energy condition is violated that 
is when (ip 2 — V)t 2 is negative. — V)t 2 is plotted for the corresponding 
numerical solutions described by fig. 9 and fig. 10. 

Figure 12: 

Evolution of a Bianchi I solution with A > y/E. 
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Fig.6 
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Fig. 10 
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